Prospects of inflation in delicate D-brane cosmology 
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We study D-brane inflation in a warped conifold background that includes brane-position depen- 
dent corrections for the nonperturbative superpotential. Instead of stabilizing the volume modulus 
X at instantaneous minima of the potential and studying the inflation dynamics with an efl^ective 
single field (radial distance between a brane and an anti-brane) (j), we investigate the multi-field 
inflation scenario involving these two fields. The two-field dynamics with the potential V{4>, x) in 
this model is significantly different from the effective single-field description in terms of the field 4> 
when the field x is integrated out. The latter picture underestimates the total number of e-foldings 
even by one order of magnitude. We show that a correct single-field description is provided by a 
field t/' obtained from a rotation in the two-field space along the background trajectory. This model 
can give a large number of e-foldings required to solve flatness and horizon problems at the expense 
of fine-tunings of model parameters. We also estimate the spectra of density perturbations and 
show that the slow-roll parameter r)^^ = M'^^V^^^ /V in terms of the rotated field t/j determines the 
spectral index of scalar metric perturbations. We find that it is generally difficult to satisfy, simul- 
taneously, both constraints of the spectral index and the COBE normalization, while the tensor to 
scalar ratio is sufficiently small to match with observations. 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 

Modern Cosmology witnessed a revolution in 1980 with 
the advent of cosmological inflation The paradigm 
has stood the test of theoretical and observational chal- 
lenges in the past two decades 0, 0] . In spite of its cos- 
mological successes to solve horizon and flatness prob- 
lems, it still remains a paradigm in search of a viable 
theoretical model. It is, therefore, not surprising that 
efforts are underway to derive inflationary models from 
string theory, a consistent quantum fleld theory around 
the Planck's scale. The discovery of nonperturbative so- 
lutions, in string theory, called D-branes y has given rise 
to new hopes in this direction. 

The D-brane cosmology is a subject of intense study 
at present. It ranges from inflation on a non-BPS brane 
because of tachyon condensation to inflation due to the 
motion of a D3-brane towards an anti-D3-brane 0, 0, 
0, [l]. Efforts also have been made to study inflation 
due to geometric tachyon arising from the motion of a 
probe brane in the background of a stack of either NS5- 
branes or the dual D5-branes . These models are 
constructed in the framework of effective field theory and 
assume an underlying mechanism for the stabilization of 
various moduli fields. 

An important step towards a realistic inflationary 
model in string theory emerged from the realization that 
background fluxes can stabilize most of the modulus 
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fields. As shown in Ref. [T3|, the fluxes in a warped 
compactification, using a Klebanov-Strassler (KS) throat 
[m , can stabilize the dilaton and complex structure mod- 
uli of type IIB string theory compactified on an oricn- 
tifold of a Calabi-Yau threefold. In fact it was demon- 
strated in Ref. [l^ that all the closed string moduli 
can be stabilized by a combination of fluxes and non- 
perturbative effects. The nonperturbative effects are 
mainly responsible for stabilizing the Kahler moduli; they 
arise, via gauge dynamics of either an Euclidean D3- 
brane or from a stack of n D7-branes wrapping super- 
symmetrically a four cycle in the warped throat. The 
warped volume of the four cycle controls the magnitude 
of the nonperturbative effect since it affects the gauge 
coupling on the D7-branes wrapping this four cycle. 

The above formalism could lead to the construction 
of a realistic inflationary model [Tsj . in string theory, 
which is built upon the compactification data (see also 
Refs. The inflaton potential can be obtained by 

performing string theoretic computations involving the 
details of the compactiflcation scheme. In this setup in- 
flation is realized by the motion of a D3-brane towards 
a distant static anti-D3-brane, placed at the tip of the 
throat. The position of the moving brane in the com- 
pactiflcation manifold is identified with the inflaton field. 
To be more precise, the location of the mobile brane can 
be labeled with five angular coordinates and one radial 
coordinate, r. The canonical inflaton field 4> is expressed 
by a constant re-scaling of this radial coordinate. 

In Ref. [iBl, the embedding of the D7-branes as given 
in (l6j was considered. It was assumed that at least one 
of the four-cycles carrying the nonperturbative effects de- 
scend down a flnite distance into the warped throat. It 
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was then shown that the presence of a D3-brane gives rise 
to a perturbation to the warp factor affecting a correction 
to the warped four cycle volume. Moreover, this correc- 
tion depends on the position of the D3-brane and thus 
the superpotential for the nonperturbative effect gets cor- 
rected by an overall position-dependent factor. The to- 
tal potential, that the mobile D3-brane experiences, is 
the sum of the potential (F-term) coming from the su- 
perpotential and the usual D-term potential contributed 
by the interaction between the D3-brane and the anti- 
D3-brane. When the corrections to the nonperturbative 
superpotential is taken into account, the issue of volume 
modulus stabilization needs to be re-analyzed. This has 
been carried out in Refs. [13, HI, [3 ^^'^ viability of 
inflation was investigated in this modified scenario. 

The model in Ref. [13, 13 described by a two-field 
potential V{(j), a) in terms of the inflaton (f) and the vol- 
ume modulus cr. If the mass of the modulus is much 
larger than the Hubble rate, the field a approximately 
evolves along the instantaneous minima determined by 
the condition = 0. One can obtain an effective single- 
field potential with respect to (j) by integrating out a in 
this way. However, we need to be careful for the fact 
that the actual trajectory is determined by the direction 
along the velocity of the fields. We shall introduce a new 
rotated field ip along the background trajectory and show 
that the usual single-field description holds for -0 and not 
for (j>. As a result the single-field description in terms of 
(j) underestimates the total amount of e-foldings even by 
one order of magnitude. This reflects that fact that the 
amount of inflation is sensitive to the slight change of the 
potential and also of model parameters. 

We shall also calculate the spectra of density per- 
turbations generated in this model. Again the spec- 
tral index n-jz is determined by the slow-roll parameter 
77^^ = M^jV^^/V instead of 77^^ = M^^V^^/V. Thus we 
show that the correct single-field description in terms of 
tp is crucially important to study the perturbation spec- 
tra as well as the background dynamics. 



II. THE D3-BRANE POTENTIAL 

In this section we review the derivation of the scalar 
potential [l3, [3 '^^ ^ mobile D3-brane. The fiuxes for 
the compactification of type IIB string theory on an ori- 
entifold of Calabi-Yau theory are chosen such that the 
internal space has a warped throat region. For example, 
the local geometry is taken to be the warped deformed 
conifold which is a subspace of complex dimension three 



in a four dimensional complex space defined by the con- 
straint X]^=i ~ ^ where zi are coordinates in the four 
dimensional complex space and e is real and corresponds 
to the deformation parameter. When the region of rele- 
vance for the D-brane inflation is chosen to lie far from 
the tip of the throat, the deformed parameter can be ne- 
glected. We can then choose Za = (21, ^2, ^3) as the three 
independent complex coordinates for the position of the 
D3-brane (open string moduli) and use the conifold con- 
straint to express 2:4 in terms of these three coordinates. 
The throat can be glued into a compact space, which is 
assumed to have a single Kahler modulus p. 
One can take the Kahler potential to be 
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so that the Kahler metric on the conifold, fc^^ ^ = dadpk^ 
is Ricci-flat. Thus the metric of the deformed conifold 
can be written as dsg — df^ -I- f^ds|,i,i, where ds^i 1 
is the metric on the Einstein space T^'^, which can be 
expressed in terms of five angular coordinates. These 
angular coordinates and the real radial coordinate r are 
the basis for the complex coordinates Zi , see Ref. (3 for 
details. 

The A^ = 1 supergravity F-term scalar potential in- 
volving the DeWolfe-Giddings Kahler potential K. and 
the super potential W is given by 




2,K^WW 



(2) 



where D^W^ = d^W 
{p, Zq,; q; = 1, 2, 3} and : 



M 



pi 



8nG. The total 



Kahler potential K. depends upon p and the open string 
moduli Za, and is given by [20| 



IC{p,p,Za,Za) = -3M^illl[p + p~jk{Za,Za)] 

= -3Mp2 mc/. 



(3) 



where k{za,Za) is the little Kahler potential, as defined 
in Eq. ^ , for the metric on the conifold. The parameter 
7 is proportional to the ratio of warped volumes of four 
cycle (before the D3-brane enters the throat) and the 
three-fold. In the expression ^ for the scalar potential, 
K,^^ is the inverse Kahler metric which can be derived 
from the Kahler potential K,. This leads to the following 
result 
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{p + P + 7(fc,a k'^^k,^ -k)} \W,p p - i{WW,p +c.c) 



+ W,pW,o. +c.c) + W,p 



(4) 



r 



The first line in Eq. @ is the standard Kachru-Kallosh- 
Linde-Trivedi (KKLT) F-term potential [Hi, while the 
rest of terms owe their existence entirely to the depen- 
dence of the nonperturbative superpotential on the open 
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string moduli [T3 (see also Refs. |2l|). 

This expression simplifies further, when we make use 
of the conifold metric (and its inverse) computed from 
the Kahler potential given in Eq. ([T]), to 
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3(W„ +C.C.) + l{z^W,pW,a. +C.C.) + -r ( S"^ + 

2 7 \ 2 r-' 



.(5) 



r 



Note that the expression for U, now, is simply U{p, r) = 
p + p — (3/2)7r^. To obtain the full potential, wc add 
to the contribution of an anti-D3-brane at the tip of 
the conifold, including its Coulomb interaction with the 
moving D3-brane p^ : 



VD{p,r) 



Do 



U^p,r) 



1 



3Z?o 



167r2(r3r2)2 



(6) 



where Dq = 2/iq"^T3 is twice the warped D3-brane tension 
at the tip of the throat. The total potential on the mobile 
brane is thus given hy V = Vf(p, Za) + Vd{p, r). 

The superpotential W — Wo + Wnp where Wq is the 
Gukov-Vafa-Witten flux super potential [l^ and Wnp ~ 
A{za)e~^'' is the contribution from the nonperturbative 
effects. Here the factor, b = 27r/ri,, arises from gauge 
dynamics on a stack of n number of D7-branes. These 
D7-branes wrap a four-cycle in the warped throat pre- 
serving supersymmetry which is specified by the embed- 
ding equation f{za) = 0. The presence of the D3-brane 
also leads to a perturbation to the warp factor which 
results in a correction to the warped four-cycle volume. 
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This correction has been found to be D3-brane position 
dependent and, in fact, is responsible for the pre- factor 
A{za) in Wnp [2^. In Ref. [l^ this position dependence 
is found to be 



A(z„) = Aq 



/(O) 



1/n 



(7) 



where A^ is related to the threshold corrections which 
depend on complex structure moduli and for us it is just 
a constant parameter as these are already stabilized by 
the flux background. Note that the relation ([7]) can be 
derived gcnerically without specif ying the embeddings 
in the warped deformed conifold [24| . The potential, 
V{p,r, Zi), thus is a complicated function of the Kahler 
modulus, the radial and the five angular coordinates of 
the mobile D3-brane. However, the angular coordinates 
can be integrated out as mentioned below. 

Choosing the holomorphic embedding of Ref. [3, 
given by the equation f{zi) = p — zi = 0, we have 
A{zi) = Ao{l — Using this and also setting 

p = (7 + tr in Eq. ^ we find 



Vf 



K%\A\'' 



-2ba 



3t/2 



2ba + 6 + 6VKoe'""Re ( — 



ibr 



3 p{zi + zi) - 2\zi\'' 



2n 



r 



67 ' 



\Zl\ 



1 



2j.3 J 7i2|^ _ 2^|5 



(8) 



Note that this potential has only one term that depends 
on r. The potential for r is minimized when this term 
takes the smallest possible value. Since the coefficient of 
this term contains Wq which is negative, integrating out 



this field amounts to replacing e'^^'^ /A by | Similarly, 
the five angular coordinates, which describe the position 
of the D3-brane on the base of the throat, are periodic 
coordinates on a compact space. Thus, the potential in 
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these fields is either constant or has discrete minima for 
some fixed values of these five fields. Since the radial 
motion of the branc is of special interest, we can focus 
on the stable trajectories in the angular directions which 
minimize the potential. 

In the case of the embedding /(zi) = 0, the stable 
minima in angular directions occur only for trajectories 
that obey zi = -r^/^/v^ From the DBI action of 
D3-br ane, on e can read out that the canonical field is 
(j) = y^3T^/2r, which corresponds to the approximation 



(when compared to the result obtained from the Kahlcr 
potential) that p + p ^ {■y/Tz)^)^ and that the field a 
does not change much. Using this expression for zi in 
Eq. ([5]) and defining the minimal radial coordinate (po- 
sition) of the D7-brane embedding to be = 2p? i.e., 
4>1 (3/2)T3(2^2)2/3^ pg^j^ ^jjp potential, in- 
volving only two fields, the field (p and the volume mod- 
ulus a, in the following form 



y(0,a) = 



h\Ao\ 



-2ba 



26ct + 6 - 6e' 



bij 



\Ao\ gi/n(0) n 



3/2 



(9) 



where 



C7(0,a) 



2(7 - . 



= 1 



T3 



= i^n 1 - 



3/2 



27Z?n 



647r2 



(10) 

(11) 
(12) 



c = l/(6^77^3</'^). 



and 

We will use this form of the two-field potential in our 
analysis in the next section keeping in mind that Aq and 
Wo have the dimension [mass]'^ while for _Do it is [mass]*. 
It should be noted that the field a is not yet in the canon- 
ical form. From the Kahler potential ([3]), we find that the 
kinetic term, —Kppdf^pd'^p, becomes canonical (at the tip 
of the throat, which coincides with the same approxima- 
tion made for the canonical field 0) if we consider the 
field X defined by 



X 
Mpi 



IntJ . 



(13) 



In what follows we shall examine the viability of inflation 
in the frame work of two-field dynamics described by the 
potential V^cj), x)- This potential amounts to replacing a 
in Eq. (HD by cxp(y273 x/Mpi). 



the equations of motion are given by 

1 



H = - 



2Ml, 



+ 2,H(j) + = , 
X + 3i/x + Kx =0, 



together with the constraint equation 



1 



1 



1 .2 

2^ 



(14) 

(15) 
(16) 



(17) 



where a dot represents a derivative with respect to cos- 
mic time t and H = a/a is the Hubble parameter. For 
later convenience we introduce the following slow-roll pa- 
rameters 



V 



V 



•n4>4> 



, r2 "^M 



^^^Pi y 



(18) 



The squared masses of the fields (j) and x a-re defined 



by 



= Vm and 



XX 



Vxx^ respectively. Since 



the Hubble parameter approximately satisfies the rela- 
tion ~ V/M^^ during inflation, the slow-roll param- 
eters ?700 and rj^-^ are given by 



III. BACKGROUND EVOLUTION 

In this section we discuss the dynamics of inflation in- 
duced by the D3-branc potential ^ but with the canon- 
ical field X replacing the field a. In the flat Fricdmann- 
Robcrtson- Walker (FRW) metric with a scale factor a 



3i?2 



Vxx 



(19) 



Since we are considering the case in which the field (/> 
plays the role of the inflaton, we require that rj^^ is not 
much larger than unity. Meanwhile the modulus field x 
can be heavy {rj^x ^ 1) or light {rjxx ^ 1) depending 
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FIG. 1: The trajectory of the scalar fields in the {(p, a) plane 
for the model parameters n = 8, Ao = 1, tJo = bao — 10.1, 
LOF = 9.9951, Wo = 3.496 x 10""*, Do = 1.215 x 10"* and 
(j>fj, = 0.25. The curves (a) and (b) correspond to the trajec- 
tories derived by solving the background equations numeri- 
cally for the initial conditions (j>i/4'^ = 0.8 and 4'i/4>ti = 0.5, 
respectively. The initial values of the field x s^re chosen to 
satisfy the relation (^(7)1 with ctq = ijJo/b — 12.8597. We also 
plot the approximate trajectory (|20|) . 



on the model parameters. The latter corresponds to the 
situation in whieh two stages of inflation can be realized. 

When the x mass is much larger than the Hubble rate 
(%x ^ 1)' field X rapidly rolls down toward instan- 
taneous minima of the potential given by V^^ = 0. In 
this case it was shown in Ref. [13, [iHl that the evolu- 
tion of the non-canonical field a = exp ( 1^/2/3 x/Afpi) is 
approximately described by 



C3/2 



3/2' 



(20) 



We note that this relation was also derived in Ref. [l9|. 
Here ctq = 'iiMlJTz and C3/2 (1 - l/2wi.) /nw_F, 
where ujp satisfies the relation 3e'^^ | Wo|/|Ao| = 2ujf + '^- 
The condition, ^ 1, is assumed to reach the expres- 
sion (|^D|) . We also note that the approximation (PO)) is 
not accurate in the large 0/0^^ region close to 1. 

In Fig. [1] we illustrate a typical example for the tra- 
jectory of the two fields in the ((/), a) plane. The case 
(a) corresponds to the initial conditions 4>i/(f>fi = 0.8 and 
CTi satisfying Eq. (PO)) . Since the approximation (|^ is 
not accurate in the large region, the initial posi- 

tion of the field a does not exactly match with the local 
minimum of the potential {V^ = 0). The field a quickly 
oscillates around instantaneous minima of the potential 
at the initial stage. In this case the system enters the 



slow-roll inflation stage around (j)/(f>fj, ^ 0.5 after the field 
X is almost stabilized at the instantaneous minima. If 
we choose smaller initial 4>/4>fi (such as the case (b) in 
Fig. [T|) , the period of the oscillation of x becomes very 
short. Figure [1] shows that the accuracy of the approxi- 
mation (|^D|) becomes better for smaller values of 0/0^. 

In the upper panel of Fig. [5] we plot the potential 
V{(f>,x) s-s a- function of 0/0^ for the same model pa- 
rameters as in Fig. [TJ Note that the field 4> is obtained 
by projecting a two-field trajectory into the (j) direction. 
Since ry^^ is of order 10'^ in this case, the mass of the 
field X is much larger than the Hubble parameter. When 
x/Mpi > 3.141 the potential does not possess instanta- 
neous minima in the region < 0/0^ < 1, but they 
appear for xA^pi < 3.141. We require that x/Mpi is 
initially smaller than 3.160 for (j)i/4>^ = 0.8 in order to 
avoid that the field evolves toward the forbidden region 

The numerical simulation in Fig. [5] corresponds to the 
initial conditions 0i/(/)^ = 0.8 and XiM^pi = 3.1385, 
which satisfy the relation ((20)l . In this case the field 4> 
slightly evolves toward larger at the initial stage, but it 
soon begins to roll down to instantaneous minima which 
move toward smaller (f> with the decrease of x- From the 
lower panel of Fig. [2] we find that the field x does not exist 
at the instantaneous minimum initially (because of the 
breakdown of the approximation (|20p ) and then evolves 
toward the instantaneous minima with oscillations. The 
upper panel of Fig. [2] shows that the instantaneous min- 
ima in the 4> direction disappear for xZ-^'^pi < 3.130. This 
is the signal of the end of inflation. Thus the field x does 
not evolve much during inflation in the above case, as 
required from the stabilization of the volume modulus. 
However, as we will see below, it is of crucial importance 
to incorporate the dynamics of the field x- 

In Ref. [l3 the authors reduce the two-field potential 
^ to the single-field one by substituting the relation 

for m-. 



K(0) = X^(0,a,(0)). 



(21) 



In Fig. Owe plot the numerically obtained potential V as 



a function of 



to compare with ((2T|) . The potential 



in the two-field system is flatter than in the case derived 
under the single (p field approximation. This implies that 
the real two-field system chooses a trajectory which gives 
a larger amount of inflation. 

In Fig. U we show the evolution of the number of e- 
foldings In a in terms of the function of 0/0^ with the ini- 
tial conditions 0^/0,, = 0.8 and 4)i/(t)^ ~ —1.0 x 10~^°?n, 
where m = lQ^'^AIp\. In the two-fleld case inflation oc- 
curs aroimd the region 0.3 ^ 4'/4'tJ. ~ 0-5, which leads 
to the number of e-foldings In a = 67 at the end of infla- 
tion. We find that this is not sensitive to the change of 
initial conditions as long as (t>i/(j)^ is larger than 0.5. The 
change of initial velocities of scalar fields hardly affects 
the evolution of the slow-roll regime, unless we choose 
unnaturally large initial velocities. In the single (j) field 
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FIG. 2: The potential 1/(0, x) for the model parameters n = 
8, ylo = 1, uo = &O-0 = 10.1, UF = 9.9951, Wo = 3.496 x IQ-*, 
Do = 1.215 X 10"* and (|)^, = 0.25. The solid curves cor- 
respond to the one derived by solving the background equa- 
tions numerically for the initial conditions cj>i/4>^ = 0.8 and 
4>i/4'iJ. ~ —1-0 X lO^^'^m, where m = 10~^Mpi is a mass to 
normalize time t. Note that the initial conditions for the field 
X are chosen to satisfy the relation (|20p . i.e., Xi/^'^pi ~ 3.1385. 
The upper panel shows the potential in terms of the function 
of (fi/cj)^ for several fixed values of x = x/^^pi (plotted as dot- 
ted curves). The lower panel shows the potential in terms of 
the function of x/Mpi for several fixed values of (p ; 



Two field 

Single ({) field approximation 




(^/ 



(|)/<|) 



FIG. 3: The solid curve represents the potential obtained 
by numerically solving the background equations in two-field 
system for the same model parameters and initial conditions 
as in Fig. [2] The dotted curve corresponds to the potential 
(|2ip under the single field approximation. 



approximation, compared to the two-field system, we ob- 
tain a much smaller value of the number of c-foldings: 
Ina = 7. 

This discrepancy reflects the fact that the background 
trajectory along which the fields evolve is not given by 
the field (j) but by the field tp satisfying the relation 



where 



^ = (cos 9)(j) + (sin 6)x , 



tan 6* = x/4' ■ 



(22) 



(23) 



Here 9 characterizes the direction along which the scalar 
fields evolve. If the field trajectory is curved as in our case 
then we have 6^0. Since s = — (sin6')(/) -I- (cos6')x = 
from Eq. ((^51) . the fields do not move to the direction 
orthogonal to ip. 

In order to find out an effective single-field trajec- 
tory we need to obtain the potential V in terms of the 
function of ip rather than (p. If the field x is frozen at 
some particular value (x — const), we just need to de- 
rive the potential with respect to cp because ip coincides 
with (p. The hybrid inflation model with the potential 
V = 1(^2 - AP/X)^ + y^cP^x^ + ^m^f M falls into 
this category: inflation occurs along the line x = 0. How- 
ever the model ^ gives a curved background trajectory, 
which means that the single-field description in terms of 
"0 is necessary to understand the dynamics of inflation 
correctly. 

The mass squared of the field tp is given by [2^ 



V 



(cos^ d)V4,4, + (sin26l)y^x + (sin^ d)V,^^ . (24) 
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I ^ Two field 



Single ({) field approximation 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
(])/ <|,/(|,^ 

FIG. 4: The evolution of the number of e-foldings in terms of 
the function of ^/(^^ for the two-field system (solid curve) and 
for the system under the single <j) field approximation (dotted 
curve). The model parameters and initial conditions are the 
same as in Fig. (2] 
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FIG. 6: The evolution of + in terms of the function of 



with three different values of Do- Other model parame- 
ters and initial conditions are chosen as in the case of Fig. [2] 
We find that -\- are smaller than 10~* during infiation 
(0.3 < <t>/4>^. < 0.5). 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 

FIG. 5: The evolution of the slow-roll parameters 77^^ and 
7700 in terms of the function of (j>/4'fi- The model parameters 
and initial conditions are chosen as in the case of Fig. [2] The 
period of infiation is determined by the condition \ri^^ \ < 1. 

Then the slow-roll parameter, 7/^^ = Mpjl/^^/y, is 

1]^^ = (cos^ 0)ri^4, + {sin29)i],f,^ + (sin^ ^')%x ■ (25) 

In Fig.[5]we plot the evolution of the slow-roll parameters 
77^^ and r]^^ for the same model parameters and initial 



conditions as in Fig. [21 The period of the inflationary 
stage is determined by the condition \r]^ip\ < 1 instead 
of the condition \r]^^\ < 1. In fact Fig. |4] shows that 
inflation occurs for 0.3 (p/cf)^ ^ 0.5, which coincides 
with the region given by the condition l?;^^, | < 1 in Fig.O 
The slow-roll parameter 77^^ is larger than unity during 
inflation, which implies that this variable is not suitable 
to describe the dynamics of inflation. In the next section 
we show that the slow- roll parameter 77^^, instead of 77^^, 
is crucially important to estimate the spectral index of 
scalar metric perturbations. 

The scalar fields approach instantaneous minima of the 
potential so close that slow-roll parameters and be- 
come much smaller than |?7^,^| during inflation, see Fig.[Sl 
It is possible to have a larger total number of c-foldings 
than the case shown in Fig. [4] by choosing slightly differ- 
ent values of Dq. For example, when Dq = 1.218 x 10~® 
instead of Dq = 1.215 x 10~*, we find that the total 
number of e-foldings becomes Ina = 148. In this case 
the field stays at instantaneous minima of the poten- 
tial for a longer time. Then we obtain smaller values 
of £0 and , which leads to a larger amount of inflation. 
This behavior is clearly seen in the numerical simulation 
of Fig. [6l Inflation ends when the slow-roll parameter 
1 771/) 1 becomes larger than unity. Meanwhile if we choose 
Dq = 1.210 X 10~^ then the total number of e-foldings 
is found to be Ina = 43, which is not sufficient to solve 
horizon and fiatness problems. When Dq = 1.220 x 10~^ 
inflation does not end because the fields are stuck at the 
local minimum of the potential. These results show how 
the number of e-foldings is sensitive to the choice of model 
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parameters. Thus we require a severe fine-tuning as it is 
typical in the context of D-brane inflation. In the field 
equations (jl5|l and (jl6p we have dropped other kinetic 
terms coming from the Kahler potential which are sup- 
pressed by ©(^/A/pi). We have carried out numerical 
simulations by taking into account those terms and found 
that the total number of c-foldings is hardly changed (the 
difference of order 0.1). Thus the results arc insensitive 
to the inclusion of such kinetic terms. 

The dynamics of inflation is mainly affected by the 
changes of the model parameters Wq and Dq. We have 
run our numerical code for many other cases by chang- 
ing Wo and Dq and found that the inflationary dynamics 
discussed above is a rather generic feature when the field 
X is massive (%x T^^^ single field approxima- 

tion using the potential (pij) tends to underestimate the 
total amount of inflation. We have also considered the 
parameter regions in which the x mass is smaller than 
the Hubble parameter (77^^ ^1)- In such cases, we nu- 
merically find that the field x typically evolves toward 
larger values unlike the estimation given in Eq. (PO)) . In 
this case the instantaneous minima of the potential tend 
to disappear in the region < (/'/(/>p < 1, which implies a 
difficulty to obtain a large number of c-foldings. Hence in 
what follows we concentrate on the case where the mass 
of the field x is much larger than the Hubble parameter. 



IV. COSMOLOGICAL PERTURBATIONS 

In this section we discuss the spectra of density pertur- 
bations generated in the model with the potential ^ . In 
two-field models of inflation the resulting density pertur- 
bations are generally different from those in a single-flcld 
model because of the presence of isocurvature (entropy) 
perturbations (see Refs. [1^, [2^ for review). Wc de- 
note the field perturbations in </> and x ^4' ^-^d 5x- 
Along and orthogonal to the background trajectory in 
two-field space it is convenient to carry out a field rota- 
tion [ii: 



Sip = (cos 9)5(1) + (sin 9)5x , 
5s = ~{smd)5(t) + {cos6)5x ■ 



(26) 
(27) 

where 6 is defined by Eq. (|23p . Here dip and 5s correspond 
to adiabatic and entropy perturbations, respectively. 

The perturbed spacetime about the FRW background 
is described by the line element 

ds^ = -{l + 2A)dt^ + 2a{d,B - S,)dx^dt 

W [(1 - 2ip)5^j + 2d,jE + h,j] dx'dx\ (28) 

where A, B, E, Lp are scalar metric perturbations and 
hij is the tensor perturbation (we omit to write vector 
perturbations). The comoving curvature perturbation TZ 
and the isocurvature perturbation S are defined by 



H5p 



S 



Hi^5x - x5c^) 



P +X^ 



(29) 



where 5p is the total density perturbation. Using the field 
perturbations and (P7|) these are simply expressed as 



n 



ip 



s 



H5s 



(30) 



where (5-0 = 5ip + ip f/H. 

The each Fourier mode for St, 
following equations of motion [26| 



5iP^ + ZHSiP. 



and 5s satisfies the 




where fc is a comoving wavenuniber, ^ — ip ^ a?H{E — 
B/a), V,^^ is defined by (pi]) and 

Vss = {smH)V^^ - {sm2e)V^^ + (cos^e)!/^^. (33) 

This shows that in the large-scale limit (fc 0) the en- 
tropy perturbation S evolves independently of the cur- 
vature perturbation, while the curvature perturbation TZ 
is sourced by the entropy perturbation as long as the 
trajectory is curved {9 ^ 0) in the field space. 

The background trajectories we discussed in the pre- 
vious section correspond to the cases in which the field 
t\} is light and the field s is heavy relative to the Hubble 
parameter during inflation. Hence the mass term 1^^^, 
in Eq. (|3ip is smaller than the order of i?^, whereas V^ss 
is larger than . When both Y m> and V^ss are smaller 
than EP' , it was shown in Ref . [13, [lH that at the Hubble 
radius crossing (fc = aEE) the cross-correlation C-ns be- 
tween curvature and isocurvature perturbations is zero at 
lowest order in slow-roll (while it docs not vanish at the 
first order). In our case {r\,ss 3> 1) the power spectrum of 
the entropy perturbation, Vis = ^ behaves in the 

usual way ij^ss — (fc/27ra)^) for j a} '3> Vss, but after 
the mass term V,ss becomes larger than k'^jc? the en- 



tropy perturbation decreases more rapidly (T^^g oc a 



This latter stage occurs even before the Hubble radius 
crossing (i.e., V^ss > k^/a^ > H^), which generally leads 
to the weaker cross-correlation at fc = aH relative to the 
case of two light scalar fields studied in Refs. [sol. Isij. 

Wc then neglect the interacting terms on the r.h.s. of 
Eq. (|31|1 for the modes fc > aH and obtain the spectrum 
Vs^^, — (-ff»/27r)^ at the Hubble exit at lowest order in 
slow-roll of the field ip (in what follows we use the symbol 
* to represent the quantities at the Hubble radius cross- 
ing, fc = aH). Hence from Eq. (j30p the power spectrum 
of the curvature perturbation at fc = aH is given by 



(fr)' 



V 



247r2eA/4 



(34) 
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where 



We also obtain the tensor to scalar ratio 



Ml, 



V 



(35) 



Here the last approximate equahty holds under the slow- 
roll approximation. Note that we used the relation ijr^ ~ 
(2/3)61^ and ~ V/Ml^. 

One can describe the evolution of perturbations after 
the Hubble exit [k < aH) by using a transfer matrix [lOl 



1 T-jis 
Tss 



7^* 



(36) 



where Tfis characterizes the correlation between curva- 
ture and isocurvature perturbations. The dimensionless 
measure of correlation is defined by 



ns 



(37) 



which is in the range \rc\ < 1 by definition. 

If the masses of two scalar fields and x small rel- 
ative to the Hubble parameter, it was found in Refs. (3^ 
that the correlation measure jrd can be close to the or- 
der of 1 in several models of two-field inflation. This 
comes from the fact that the entropy perturbation 5s is 
not suppressed after the Hubble radius crossing. Mean- 
while if one of the fields is heavy then the amplitude of 
5s exponentially decreases {\5s\ oc a""^/^), which results 
in a very weak correlation (|rc| <C 1). Since we are con- 
sidering this latter situation for the model Q, it is a 
good approximation to neglect the correlation after the 
Hubble radius crossing. 

Using this property, the power spectrum of the curva- 
ture perturbation at the end of inflation is given by 



V 



247r2eM4 



(38) 



which holds for \rc\ <C 1. When the correlation is strong, 
the r.h.s. of Eq. is multiplied by the factor 1/(1 -r^) 
[sot . The spectrum index of the curvature perturbation, 
nji = 1 + dlnT'Tj/dlnfc, is given by 



~ 1 - 6e + 2??^^ , 



(39) 



where the slow-roll parameters are evaluated at the Hub- 
ble exit. 

The tensor perturbation hij is decoupled from the 
scalar perturbation and is frozen after the Hubble radius 
crossing. Thus its power spectrum is given by [sot 



with the spectral index 

dlnPr 



2K 



3^2M4 ' 



riT 



dlnfc 



-2e. 



(40) 



(41) 



Vt 



— ~ 16e ~ 16(60 + e^) , 



(42) 



which is valid for jrd <C 1. In the presence of the cor- 
relation the r.h.s. of Eq. (|42|) is multiplied by the factor 

As we mentioned in the previous section, the slow- 
roll parameters and for the model ([9]) are typi- 
cally very much smaller than 1 because the scalar flelds 
evolve around instantaneous minima of the potential (see 
Fig. [5]) . This implies that the tensor to scalar ratio given 
in Eq. (j42p is much smaller than one. In the numerical 
simulation of Fig. [21 for example, we obtain r of the order 
10~^ on cosmologically relevant scales. Hence this model 
satisfies the present observational upper bound: r < 0.3 
[1, [131 ■ We also note that the spectral index nx of tensor 
perturbations is very close to scale-invariant. 

Since the model (O generally satisfies the relation e <C 
|?7^^| during inflation except for the vicinity at 77^^ = 0, 
the spectral index ([39]) of scalar perturbations yields 

n-jzc^ 1 + 2r]^^ . (43) 

It is important to note that n-jz is determined by rj^^ 
instead of r/^^. This reflects the fact that the background 
trajectory is not along the (j) direction but along the ip 
direction. 

The expression P5)) is valid as long as the slow-roll pa- 
rameter 770^ is smaller than the order of unity. In the case 
of Fig. [5l for example, this can be used for 77^^ < 0.2- 
0.3 after the fields almost stop oscillating and enter the 
inflationary stage. We have also computed the correla- 
tion measure Tc numerically and have confirmed that Tc 
is very small for the modes that crossed the Hubble ra- 
dius during slow-roll inflation so that the expression (j43p 
is valid. 

In Fig. [7] we plot n-ji in terms of the function of the 
number of e-foldings N from the end of inflation for sev- 
eral different cases. For larger Dq the total number of e- 
foldings increases, which leads to the change of the curves 
from (a) to (c). The curve (b) corresponds to the case in 
which the number of e-foldings exceeds the typical COBE 
scale value iV = 60, while in the case of the curve (a) it 
does not reach cosmologically relevant scales. In the case 
(b) we obtain the value n-ji ~ 1.6 for iV = 60, which 
is too large to satisfy recent observational constraints: 
= 0.97-1.21 [H.i 

It is possible to realize the red-tilted spectrum [n-ji < 
1) that satisfies observational constraints if cosmologi- 
cally relevant scales are in the negative 77^^ region. In 
the case (c) of Fig. [71 for example, we obtain n-ji ~ 0.98 
for N = 60. It is clear from Fig. [7] that the spectral in- 
dices around iV = 60 get smaller for larger total number 
of e-foldings. 



^ Without the running of the spectral index the constraint on nji 
is severer: n-j^ = 0.93-1.01. 
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FIG. 7: The spectral index nv. as a function of the number of 
e-foldings A*' from the end of inflation. Each case corresponds 
to (a) Do = 1.213 X 10"'', (b) Do = 1.215 x 10"'', and (c) 
Other model parameters are the same as 



Do = 1.218 X 10 
in Fig. [2] 



FIG. 8: The amplitude of the power spectrum V-r. as a func- 
tion of the number of e-foldings TV from the end of infla- 
tion. Each case corresponds to (a) Do = 1.213 x 10"'', (b) 
Do = 1.215 X 10"*, and (c) Do = 1.218 x 10"". Other model 
parameters are the same as in Fig. [2] 



The model also needs to satisfy the condition of the 
COBE normalization Q 

Vn ~ 2.4 X 10"^ , (44) 

on cosmologically relevant scales observed by COBE. In 
Fig. [5] we find that in the case (b) this condition can 
be satisfied around the scale A'^ = 60, while the spectral 
index n-n is larger than observationally allowed values. 
Wc obtain the smaller n-n consistent with observations 
in the case (c), but the amplitude V-r, becomes much 
larger than the value (|^^ . The increase of the amplitude 
reflects the fact that for larger total number of e-foldings 
the fields arc stuck around instantaneous minima of the 
potential for a longer time so that e tends to be smaller 
in Eq. (j38p . Thus for the model parameters chosen in 
Figs. [7] and [8] the model does not satisfy, simultaneously, 
both observational constraints of rm and Vti- 

We have tried many other cases and have not found a 
case in which both n-ji and V-n satisfy observational con- 
straints. The behavior shown in Figs. [7] and [8] is typical 
in our model. In order to obtain a nearly scale-invariant 
spectrum, we need to adjust that cosmologically relevant 
scales (N ^ 60) exist in the region around the maximum 
of the power spectrum. In this case, however, the am- 
plitude V-n tends to be too large to satisfy the COBE 
normalization. Nevertheless we have not searched for all 
parameter spaces; presumably there may be some regions 
or isolated points in the parameter space that satisfy both 
constraints of tt-tj and V-r, simultaneously. However the 
difficulty to find such viable parameters implies that the 
model requires severe fine-tunings in order to be consis- 



tent with observations. 



V. CONCLUSIONS 

In this paper wc have studied inflation in "delicate" D- 
brane cosmology [13, [111 that takes into account effects 
of the moduli stabilization. The presence of the D3-brane 
in a warped conifold background gives rise to a correc- 
tion to the warped four cycle volume, which leads to a 
modification to the nonperturbative superpotential W . 
Hence the total potential V depends upon not only the 
inflaton field but also the (canonical) volume modulus 
field X- 

Since the field x evolves in two-field space even if the 
change is small, we need to study the dynamics of multi- 
field inflation in an appropriate way. If the mass of x is 
much larger than the Hubble parameter, this field quickly 
evolves toward instantaneous minima of the potential. 
In Ref. [13, [H, [H the approximate relation with 
respect to the non-canonical field u = exp ( ■\/2/ 3 x /-^^pi ) 
was derived under the condition = 0. This relation 
tends to be accurate for smaller <f>, as we have shown in 
Fig. [TJ The effective single field potential in terms of 
(f) can be obtained by substituting Eq. for the two- 
field potential ([9]). Numerically we have found that this 
description underestimates the total number of e-foldings 
even by one order of magnitude, see Fig. [H 

This discrepancy comes from the fact that the actual 
background trajectory is not given by the effective single 
field but by the field -0 defined by Eq. If the 
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scalar fields evolve along a straight line with constant x 
during inflation, as in the case of hybrid inflation with the 
potential V = ^{x^ - M"^ jXf + \g^^'^x^ + \m^<i?., the 
single-fleld description with respect to </> is valid because 
the V' direction coincides with the direction. If the 
trajectory is curved as in the model Q , one has to study 
the dynamics of inflation along the "0 direction. In fact in 
the numerical simulation of Fig. [5] the period of inflation 
(0.3 < 0/0^ < 0.5) coincides with the one derived by the 
slow- roll condition \'r]^^\ < 1, but the slow-roll parameter 
7^00 is larger than unity during this epoch. 

We also find that the total number of e-foldings is very 
sensitive to a slight change of model parameters. This is 
related to the fact that the naive single (j) field descrip- 
tion grossly underestimates the total amount of inflation. 
Although the change of the fleld x during inflation is typi- 
cally small, the evolution of the volume modulus nontriv- 
ially affects the dynamics of inflation. It is interesting to 
note that in two-fleld model ^ one can realize a large 
number of e-foldings to solve horizon and flatness prob- 
lems even if it is difficult in single-field context. 

We have also evaluated the power spectra of density 
perturbations generated in this model. In two-field mod- 
els of inflation there exist isocurvature perturbations in 
addition to curvature perturbations in general. When 
one of the fields is much heavier than the Hubble rate, 
the correlation between curvature and isocurvature per- 
turbations is generally very small. In this case the power 
spectrum V-jz of scalar metric perturbation is estimated 
by Eq. p8|) together with the spectral index n-fi given 
by Eq. P5|) . It is important to recognize that the spec- 
tral index is determined by the slow-roll parameter 77^^ 
instead of rj^^. We have found that it is generally dif- 
ficult to satisfy observational constraints of both the 
spectral index and the COBE normalization simultane- 



ously. This comes from the fact that as the spectrum 
approaches scale- invariant (ji-n = 1) the amplitude V-r, 
tends to be larger than the COBE normalized value 
{Viz ~ 2.4 X 10~^) on cosmologically relevant scales, see 
Figs. [7] and m 

The model we have studied is in fact "delicate" and 
needs severe fine-tunings of model parameters to satisfy 
several constraints discussed in this paper. The imple- 
mentation of corrections due to the D3-brane motion to 
multi-throat D-brane models [s^ can possibly improve 
the situation. It is to be noted that a viable model of in- 
fiation should be followed by a successful reheating [s^l 
which is a challenging problem for D-brane cosmology 
at present. The multi-throat models have recently given 
some hopes in this direction, see Refs. [s^ on the related 
theme. These are important issues of string cosmology 
which require further investigation. 
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